Spontaneous Symmetry Breaking of Population between Two Dynamic Attractors in 
a Driven Atomic Trap: Ising-class Phase Transition 
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We have observed spontaneous symmetry breaking of atomic populations in the dynamic phase- 
space double-potential system, which is produced in the parametrically driven magneto-optical trap 
of atoms. We find that the system exhibits similar characteristics of the Ising-class phase transition 
and the critical value of the control parameter, which is the total atomic number, can be calculated. 
In particular, the collective effect of the laser shadow becomes dominant at large atomic number, 
which is responsible for the population asymmetry of the dynamic two-state system. This study 
may be useful for investigation of dynamic phase transition and temporal behaviour of critical 
phenomena. 
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The phenomena of symmetry breaking, widespread in 
nature with examples from cosmology to biology, have 
been much studied 0, S 0- Recently in a vibro- 
fiuidized granular gas, spontaneous symmetry breaking 
(SSB) of temperature and population between two com- 
partments connected by a hole was reported and under- 
stood in terms of the density-dependent inelastic collision 
rates 4, 5]. Moreover there have been many works on 
fluctuation-induced transitions in equilibrium 6, 7, 8] as 
well as far from equilibrium |3, [lO: HH • The double well 
structure of these systems is very similar to the two com- 
partments of the granular box, where SSB was observed. 
In particular, we have recently studied the atomic pop- 
ulation transition between two dynamic phase-space at- 
tractors available in the parametrically driven magneto- 
optical trap (MOT) system 

In this Letter, we report on experimental as well as 
theoretical investigation of SSB of the atomic popula- 
tion between two dynamic states of the driven MOT. 
We have found that the control parameter for SSB is 
the total number of atoms in both states: The popula- 
tion equality between the two equivalent states is broken 
spontaneously above a critical number of atoms. This 
phenomenon can be well understood as the Ising-class 
phase transition. We have measured the critical num- 
ber under various experimental parameters and also ob- 
served the temporal evolution from symmetric to asym- 
metric states above threshold. The SSB mechanism is 
described qualitatively by considering two collective in- 
teractions occurring at large atomic number ^he shadow 
effect and the reradiation effect (13, ll^? llJ? H^j • Iii Par- 
ticular, the measured critical numbers are in good agree- 
ment with the analytical and the simulational results. 

The experimental scheme is similar to those reported 
in previous works on parametrically modulated MOT 
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pH -> where we observed parametric excitation, limit 
cycles (dynamic phase-space attractors), super-critical 
and sub-critical bifurcation. In particular, the bifurca- 
tions were explained by atomic double- and triple-well 
potentials in the rotating phase space. Due to fluctu- 
ating atomic motions resulting from spontaneous emis- 
sions, population transfer occurs between the two states 
of dynamic double well, which tends to equalize the pop- 
ulation of each state [Fig. ^a)]. This atomic transition 
between the two states oscillating in position space (Fig. 
1) was confirmed by observing the temporal recovery of 
the population symmetry after emptying one state. In 
this case, the recovering rates are equivalent to the tran- 
sition rates 11]. 

It is interesting to observe that the population symme- 
try, which is equivalent to zero spontaneous magnetiza- 
tion in the Ising spin system, is only maintained below a 
certain critical value of the total atomic number. Above 
the critical number, however, we have observed SSB of 
atomic population, as shown in Fig. ^b). The SSB can 
be observed under wide experimental conditions of modu- 
lation frequency / and amplitude /i, from super-critical to 
sub-critical bifurcation regions. The atomic populations 
were simultaneously measured by resonant absorption of 
a weak probe laser. The typical experimental parameters 
are as follows: magnetic-field gradient along the atomic 
oscillation direction (z-direction) b = 14 G/cm, cooling 
laser detuning J = -2.6 F, and laser intensity in the z-axis 
Iz = 0.039 Is {Is is the averaged saturation intensity, 3.78 
mW/cm^). The intensity on the transverse axes is typi- 
cally 5 times larger than that of the z-axis. The measured 
trap frequency is 43.6 (± 2.4) Hz whereas the damping 
coefficient is 160.4 (± 33) s~^, which is about three times 
larger than that expected in the Doppler theory 17]. 

Figure El^a) presents the normalized population dif- 
ference between the two dynamic states (1 and 2 with 
Ni > N2), Ap = {Ni - N2)/Nt versus the total atomic 
number, Nt = Ni -\- N2. As shown in the figure, the 
main control parameter of SSB is Nt so that SSB (or 
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FIG. 1: Snapshot images of atoms in two dynamic attrac- 
tors (a) before SSB of atomic population and (b) after SSB, 
taken by a charge- coupled- device (CCD). The total number 
of atoms is (a) 6.1x10^ and (b) 6.9x10^, respectively. The 
relative population difference in (b) is 0.63. Here / = 96 Hz, 
h = 0.9, and the abscissas are in unit of mm. 



the Ising-class phase transition) occurs above the critical 
number Nc. We have measured Ni and 7V2 by two inde- 
pendent methods: CCD images (filled black boxes) and 
probe absorption (empty boxes). Nt was varied by ad- 
justing the intensity of the repumping laser while all the 
other trap parameters remain fixed. Note that we did 
not find any other control parameters other than Nt: 
for instance, the intensity imbalance between the -\-z and 
—z laser beams did not contribute to SSB for the imbal- 
ance of up to 20%, beyond which the atomic limit cycle 
motions were not sustained. 

We have measured the critical number Nc at various 
experimental parameters of / and h. For example, at 
h = 0.86, Nc decreases gently from 7.9 x lO'^ to 4.1 x lO'^ 
as / increases from 1.95/o to 2.4/o (/o is the MOT trap 
frequency along the z-axis). At / = 2.1/o (= 90 Hz), 
Nc also decreases from 11.9 x 10^ to 6.2 x 10^ as h in- 
creases from 0.64 to 0.86. In brief, Nc increases with 
the transition rate W: when f or h increases, W de- 
creases (see Ref. ^ for details), and consequently Nc 
becomes decreased. However, SSB is not observed out- 
side the above region of parameters, that is, near the 
super-critical or sub-critical bifurcation points. Around 
the super-critical bifurcation point, W becomes too large 
to load atoms enough to produce SSB in our experimen- 
tal system. Near the sub-critical bifurcation point, on 
the other hand, despite the low Nc^ only the population 
of the central stationary state among the triple wells in- 
creases, whereas those of the two dynamic states do not 
mcrease above Nc fll^. 

In order to have a better understanding of SSB, we 
have investigated the temporal evolution of the popula- 
tions of each state when SSB occurs. Figure [2fb) shows 
the atomic populations recorded by the absorption of 
the probe laser. As can be found, in the initial load- 
ing stage, the number of atoms in each state increases 
at the same rate and their growth is indistinguishable 
with each other. As the loading process is finished at 
about 20 s elapse, however, the population of one state 
increases whereas the other state is depopulated. The 
fact that the total atomic number is conserved within ex- 
perimental errors during the SSB process indicates that 
SSB originates not from different loading rates to each 
state but from the transfer of atoms from one state to 
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FIG. 2: (a) Experimental data of the normalized popula- 
tion difference Ap versus the total atomic number Nt- Here 
/ = 88 Hz (= 2.0/0) and h = 0.86. The solid curve is the fit- 
ting by the Ising model function with Nc — 6.1 (±0.4) x 10^. 
(b) Temporal evolution of SSB. The split theoretical curves 
represent the normalized populations of the two dynamic 
states [Nt = 6.5 x lO'^ and Ap = 0.52). The simple sum 
of each population, i.e., the normalized total atomic number, 
is also shown on the top. 



the other. Moreover, when we place a kicking laser near 
the center of the two dynamic states in order to block any 
transitions between the states, the population symmetry 
is recovered. These evidences confirm that SSB occurs 
due to the atomic transfer between the two states. 

Based on the fact that SSB appears above the critical 
number, one may conjecture that the underlying mecha- 
nism of SSB is related to the collective effects of atoms 
occurring between the two dynamic states ; 1 3 . 11^ 113 . IT5| . 
There are two such collective mechanisms associated with 
the MOT atoms. One is the shadow effect caused by ab- 
sorption of the cooling lasers in the z-axis due to atoms 
in one of the two states, which results in the reduction 
of the laser intensity for atoms in the other state. The 
other is the reradiation effect that arises when an atom 
reabsorbs photons that are spontaneously emitted by an- 
other atom, which produces the repulsive Coulomb-like 
forces between the two atoms. 

The reradiation effect, in fact, contributes as an ob- 
stacle to SSB. As the number of atoms in one of the 
two states becomes dominated due to fluctuations, the 
repulsive reradiation force becomes bigger outside the 
more-populated atomic cloud. As a result, this effect 
prevents atoms in the smaller-number state from being 
transferred to the larger-number state, which results in 
the recovery of population symmetry between the two 
states. On the other hand, the shadow effect acceler- 
ates atoms to move from the smaller-number state to the 
larger- number state: due to the bigger shadow effects 
associated with the larger-number state, the net atomic 
force is directed toward the larger- number state, which 
enhances SSB further. 

For a theoretical understanding of SSB process, we 
have adopted the phase-space Hamiltonian-function for- 
malism developed in Ref. ^ to account for the tran- 
sitions between the dynamic double wells [ll|. We then 
have generalized the approach to include the shadow ef- 
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feet as well as the reradiation force. This approach pro- 
vides quantitative analysis of nearly all the fundamental 
characteristics of SSB such as the critical number and 
the temporal evolution. 
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FIG. 3: (a) Hamiltonian function H^{Xi,Yi) [Eq. (0] for 
the symmetric state without collective effects, (b) Hi + iff 
with the shadow effect included, which leads to SSB in the 
phase-space potential, (c) Reradiation interaction if/^, which 
opposes SSB. (d) H^ + Hf + H^, which shows slightly reduced 
SSB with respect to (b). Here iVi = 3 x 10^ A^2 = 1 x 10^ 
I, = 0.035/s, It = 8h,b= 14 G/cm, 6 = -2.5r, / = 2/o, 
h = 0.9, and the transverse (longitudinal) spatial width of the 
atomic cloud is Rp — 1 mm (Rz — 2 mm). 

From the Doppler equation of MOT pj] . IT^ , one can 

derive the phase-space Hamiltonian function of an i-th 
atom without the collective effects as, 

H^{Xi, Yi) = \{^- l)X^ + \{^ + l)Yl - + Y^ f 

where fi = 2{f — 2/o)//i/o, Xi and Yi are the two scaled 
canonical variables in the rotating phase space, as repre- 
sented in Fig. [3Ja). The total Hamiltonian Hq without 
the interaction terms is just the summation of each Hf, 
Ho = Y^fj^i H^' Note that, in our analysis, the two os- 
cillating atomic states can be approximated as a static 
double well in the phase space [Fig. Ol^a)]. For example, 
the right- and left-side cloud with respect to the center 
of the limit cycle motion along the z-axis corresponds 
to the positive and negative state, respectively, in the Y 
axis at a given modulation phase of 0°. Note also that 
the maximum points in the phase-space potential in Fig. 
01 indicate attractors. 

Let us first consider the shadow effect that is respon- 
sible for SSB. At the specific modulation phase of 0°, 
we assume the right-side atomic cloud is the state 1 and 
the left-side cloud is the state 2. Because of the Zeeman 
shift, atoms in state 2 and state 1 absorb preferentially 
the cooling lasers propagating in the +z and —z direc- 
tion, respectively. Now we consider another j-th atom in 
state 1 or 2, which absorbs the laser by the amount = 
(j^j^UpIz^ where a;^ is the absorption cross-section of the 



j-th atom and Up is the density of atoms in the x?/-plane. 
One can then easily find that each laser photon absorbed 
by the j-th atom effectively results in a cooperative force 
(or acceleration) on the i-ih atom concerned, whose mag- 
nitude is given by Cscr]^nplz/ Is- That is, when the atom 
j is in state 1 (i.e., in the right-side cloud absorbing the 
photons propagating in the -z axis), the direction of the 
effective force experienced by the i-th atom is positive 
along the z-axis, whereas it is negative when the j-th 
atom is in state 2. If one considers, for convenience, the 
i-th atom is near the center, the net effective force exerted 
on the i-th atom is given by Csl\/Is = (^i — ^2) 
C si a/ Is- Therefore the i-th atom is transferred to the 
larger-number state on the right (state 1) [Fig. O^b)]. 
Note that if one considers the tt phase of modulation, 
although the location of state 1 (2) is now exchanged 
to the left (right), the net force is still directed to the 
larger- number state of 1. 

The Hamiltonian Hf for the shadow effect is then de- 
rived as, when summed over j in the states 1 and 2, 

iei,2 

= a{N^-N2)Yi + {N^+N2)H' , (2) 

where Uj = a = 2Cs<yLlz/Is''^^PC^^^vf^''^Iil (<^l as- 
sumed independent of j), C's = TikT/2m{l + 4(5^/r^), 
C = 27rhf§/pf, r, = ^A7r(3ff§/3Ao{(3^ + iTr^f§), P is 
the damping coefficient, and Aq is the coefficient of the 
third-order term in the Doppler equation of MOT 1^. 
H^ is a given coefficient that is practically independent 
of j, with no contribution to SSB. Here (Jl is regarded as 
having no dependence on velocity and position of atoms, 
which are assumed uniformly distributed in the x?/-plane. 
We will discuss later about more realistic treatment of the 
shadow effect with Monte- Carlo simulations. As shown 
in Fig. Ol^b), the shadow effect makes the potential of 
the larger- number state (state 1) deeper, whereas that of 
the smaller- number state (state 2) shallower. As a result, 
more atoms will be transferred from state 2 to state 1, re- 
sulting in SSB of the atomic population. In fact, however, 
there are competitions between the shadow-induced SSB 
and the fluctuation-induced symmetry-recovering transi- 
tion. Therefore the critical number is determined by the 
balance between the shadow effect and the diffusion. 

Let us now consider the symmetry-preserving reradia- 
tion effect. Figure EJc) presents the results of reradiation 
interaction i^/^, which reduces the SSB effect (detailed 
expression of H[^ will be given elsewhere). Briefly speak- 
ing, the reradiation effect increases the critical number by 
reducing a to a — Cr, where Cr = /Tcr|/47rc[37r + 8(4 + 
31n2)]/67r^i?2/^C^^^/- The calculation also shows that, 
if the transverse laser intensity is over 10 times larger 
than that of the z-direction laser, the reradiation effect 
dominates over the shadow effect, which inhibits SSB for 
every Nt. In practice, we have experimentally observed 
that the recovery of symmetry appears at about 20 times 
the z-laser intensity, which is a strong and independent 
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evidence that the reradiation hinders SSB. 
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FIG. 4: (a) Ap vs Nt, obtained by Monte- Carlo simulations 
with 10^ atoms. The solid curve is a plot of Eq. @- (b) 
Simulation curves for the temporal evolution of SSB, in good 
agreement with the experimental data (Fig. |2Ib)). Here Wo 
= 1.0 s-^ and Nt = l.lA^c 

Let us discuss the temporal evolution of SSB, which 
can be described by the simple rate equation, dAp/dt = 
Wi2{l + Ap) + W2i(l - Ap), where the transition 
rate Wi2(2i) from state 1 (2) to state 2 (1) is 
WoeMTApNr/Nc) and = D/2aC ^JTTTfifi). 
Here Wo is the atomic transition rate without the col- 
lective effects, D is the phase-space diffusion constant 
[Tl| . and /(/i) is an 0(1) function 18]. Interestingly, the 
above rate equation leads to the steady-state solution 
given by 

Ap = tdinh{ApNT/Nc) . (3) 

This is a representative equation of Ising-class phase 
transition, which is plotted in Figs. [2fa) and 4(a). 

To manifest further the relation with the Ising model, 
let us consider a simplified model where each atom i 
either belongs to state 1 {Yi = ^/i + 1) or to state 2 
{Yi = —\/fi + 1). The activation energy Si due to the 



shadow effect of the j-th atom is —2(f{fi)ajYi. The 
total interaction energy is thus X^iXlj^^^i = "(^/2) 
(A^i - 7V2)^ where J = 2aU{lj)^fi^^' The free en- 
ergy of this model system is then F = —{J/2){NtAp)'^ 
+ DNt{[{1 + Ap)/2] ln[(l + Ap)/2] + [(1 - Ap)/2] ln[(l - 
Ap)/2]}. The equilibrium value of Ap is determined 
by the condition dF/dAp — 0, which results in Ap = 
tanh(ApA^T-^/^) that is exactly the same as Eq. (|2| 
with Nc = D/J. The simple analytical values of Nc are 
in qualitative agreement with the experimental results of 
Fig. Ha). 

We also have performed Monte-Carlo simulations with 
more realistic consideration of the shadow effect and the 
reradiation force: we included the dependence of (Jl on 
the position and velocity, the transverse laser-intensity 
profile, and the random forces due to spontaneous emis- 
sions. Figure El a) shows Ap versus Nt^ which is very 
similar to the experimental results in Fig. Efa). Figure 
^b) presents the simulation curves for the temporal evo- 
lution of SSB. Here we have just included the shadow 
effect and used a diffusion constant that is 2.5 times the 
value derived from the simple Doppler theory. When the 
reradiation force is included in the simulations, however, 
Nc is slightly increased and SSB does not occur if the 
transverse laser intensity is over 10 times the z-laser in- 
tensity. In conclusion, nonlinear dynamic study of driven 
cold atoms may be useful for dynamic phase transition 
and temporal dependence of critical phenomena. 
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